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Abstract. These lectures include an introduction to the partonic description of the 
proton, the photon and the 'colour singlet', as seen in inclusive and semi-inclusive DIS, 
in e + e~ collisions, and in diffractive processes, respectively. Their formal treatment 
^vq , using structure, fragmentation, and fracture functions is outlined giving an insight into 

' the perturbative QCD framework for these functions. Examples and comparisons with 

experimental data from LEP, HERA, and Tevatron are also covered. 

os 

On 

p,; INTRODUCTION 

The discovery of asymptotic freedom, one of the most significant properties of 
strong interactions embodied by Quantum Chromodynamics (QCD) opened, more 
than 25 years ago, a new chapter in our understanding of the structure of matter 
^ ■ which has has been actively followed by theoreticians and experimentalists ever 

since. The short distance structure of hadrons, together with the production of jets 
in hadronic collisions are paradigmatic among the strong interaction phenomena 
successfully accounted by QCD and even though the standing of the theory is today 
well established, further theoretical refinements and the corresponding experimental 
validation renew constantly the original enthusiasm of the high energy physics 
community. 

These lectures intend to provide an overview of the more recent topics of high 
energy collisions related in a way or another to perturbative QCD. First we briefly 
remind the essentials of QCD, including the main features of partons (the quarks 
and gluons), which are the true protagonists in the story. Then, we refer to what 
is known about the partonic structure of three of the main benchmarks of QCD, 
the proton, the photon, and the singlet colour or 'pomeron'. Finally we will try to 
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draw the connections between their corresponding structures which in some way 
relate the physics made in the three main HEP laboratories. 

As usual, many interesting and highly active topics have been excluded from 
the lectures in favor of a more detailed analysis of the covered points. These 
include, for example, those related to the spin structure of the proton, which have 
driven an ongoing series of polarized experiments and a great deal of theoretical 
discussions [1]; heavy flavours, which involve very subtle theoretical approaches, 
and perturbative QCD beyond NLO, which is relevant for the most recent high 
precision experiments [2]. 

QCD 

The strong interactions among quarks and gluons are described by Quantum 
Chromodynamics (QCD), the non-abelian gauge theory based on the gauge group 
SU(3)c- Each quark flavour corresponds to a colour triplet in the fundamental 
representation of SU(3) and the gauge fields needed to maintain the gauge symme- 
try, the gluons, are in the adjoint representation of dimension 8. Gauge invariance 
ensures that gluons are massless. The QCD Lagrangian may be written as 

£qcd = -\ F; u Fr + i>i (< r D, - m) ^ (1) 

where 

F a , u = d,G a u - d u Gl + g f abc G b , G cv (2) 

stands for the gluon field tensor, are the quark fields and the covariant derivative 
is defined by 

D, = d,-igT a Gl 

The strong coupling is represented by g and indices are summed over a = 1 , . . . , 8 
and over i = 1,2,3. Finally, T a = A a /2 and f a b c are the SU(3) generators and 
structure constants, respectively, which are related by [T a , T&] = i f a b c T c . 

Like in Quantum Electrodynamics (QED), the procedure employed to deal con- 
sistently with the divergences that occur in the computation of strong interactions 
beyond the tree level, shows that the actual strength of the QCD coupling depends 
on the energy scale of the process. But in opposition to QED, this renormalized 
strong coupling is small at high energy (momentum) , going to zero logarithmically, 
i.e. QCD has the property of asymptotic freedom. Consequently, in this regime 
perturbation theory is valid and tests against experimental data can be performed 
in terms of hadrons. Figure 1 summarize the basic QCD perturbative processes 
appearing in different circumstances. 

Experiments with e~ e + colliders provide clean results for QCD tests. Recently, 
a huge amount of experimental data came from the HERA electron-proton collider 
[3] and also from the Tevatron at Fermilab [4]. In both cases, there is a hadronic 
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FIGURE 1. Basic processes in perturbative QCD. 

remnant that make the analysis a little more involved. All this experimental evi- 
dence support the existence of quarks being colour triplets of spin 1/2 and of gluons 
being vector octets. Moreover, the presence of the QCD coupling has manifested 
itself in different measurements, as well as the above mentioned property of asymp- 
totic freedom. This information comes mainly from the study of the so called two- 
and three-jets events [5]. 

When a given process needs a higher order in perturbation theory to be known, 
it is necessary to compute not only the renormalized strong coupling constant 
but also the appropriate corrections to the relevant cross-sections. As is usual in 
Quantum Field Theory, a regularization-renormalization procedure is in order, just 
to absorb divergences into the definition of physical quantities. This prescription 
requires the introduction of a new scale //, fixing the renormalization point, and all 
renormalized quantities begin to depend on it. Nevertheless, different prescriptions 
must end with the same predictions for observables. 

In order to illustrate how the general procedure works, ending with the Renor- 
malization Group equations that guarantee that physical observables do not depend 
on the scale /i, let us show what happens with Green functions. Just to remember 
the procedure, let us begin with a single particle irreducible Green function T. In 
general, to control divergences, one has to introduce an ultra-violet cut-off A, or 
the equivalent dimensional regularization parameter, in the loop momentum inte- 
gral defining the T. In a renormalizable theory, as QCD is, a renormalized Green 
function is defined as 

rit(pi,g,fj) = Z T (g ,A./fi) Tu(p h g^K) 

where Pi stands for the external particle momenta, go and g are the bare and the 
renormalized couplings, respectively. This T R is then finite in the limit A — > oo but 
it depends on the scale at which the value of the renormalized quantities are fixed, 
the prescription parameter /x. The function Z-p is a product of renormalization 
factors. Due to the fact that the unrenormalized IV is obviously independent of /x, 
one has to demand 

a/x 



and consequently, the Renormalization Group Equation (RGE) 

( /i |l + /9 ^ + 7 ) r «(^,/i) = (3) 

has to be verified. Here 7 is the anomalous dimension, depending on the particular 
Green function under consideration, and the beta-function is universal 

P dZ r dg 

1 = Y V ^ ^ ) = /i ^ (4) 

If there is only one large momentum scale Q, or Q 2 as it is standard to quote, as 
it is the case here, one can express all pi in terms of a fixed fraction Xi of Q. Then, 
defining the so called evolution variable 



it is possible to introduce the momentum dependent, or running coupling through 
the integral 

rff(*) dg' 



9{t) dg'J^l 

9(0) P(g') 



I g(o) P(g') 

and the general solution of the RGE reads 

T{t,g{0),Xi) = T{0,g{t), Xi ) exp 



This solution explicitly shows that the Q-scale dependence of T arises entirely 
through the running coupling g(t). Introducing now the usual notation 



,2 



Qs = ^ (6) 



one can expand the beta-function in a power series in a s 

da s _ po 2 f3i 3 



# a ') = ^ = -^ a '- 4^2 (7) 
where it results [6] that 
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fa = ll- -N f /h = 51-—N f (8) 

Here Nf indicates the number of flavours that can be excited (with mass less than 
fi) at the scale //. 

It is clear that the solution of the differential equation for a s introduces a con- 
stant, called Aqcd, which has to be fixed by using experimental data. The resulting 
a s can be written as 




FIGURE 2. QCD running coupling. 



, a v 4vr f 2 A ln[ln(/i 2 /A^)] ) 

Ago*) - Pq Hfi y AhcD) I 1 " ^ ln(^/A^) J (9) 

This expression for the running coupling shows clearly the property of asymptotic 
freedom of QCD, i.e., the coupling vanishes when the scale becomes asymptotic, 
namely /i — > oo. Consequently, in this momentum regime, perturbation theory is 
valid. 

A very clear quantitative test of perturbative QCD is provided by the measure- 
ment of a s in different processes at different scales Q 2 . In Figure 2 there is a 
summary of the various determinations of a s [7]. 

The present world average [8] for the coupling at the Z° mass is 

a s (M z ) =0.119 ±0.002 

which implies 

Aqcd = 220 + 78 - 63 MeV 



corresponding to five flavours excited and in the conventional MS prescription 
commonly used. 



PROTON 

Having reviewed the essentials of QCD, we now proceed with the analysis of 
one of the most rich benchmarks of perturbative QCD which is the short distance 
structure of the proton. 

We start this section giving a set of definitions of the commonly used relativistic 
invariants related to Deep Inelastic Scattering, DIS, one of the fundamental ex- 
perimental tools for hadron analysis, together with the corresponding formulae for 



the neutral and charged current cross-sections, where the structure functions are 
introduced. These functions will be latter expressed in terms of the quark-parton 
model including QCD corrections. For a more detailed treatment see [5] 



Lorentz Invariants 



The scattering of a lepton from a proton (or in general a nucleon, a hadron 
or a nucleus), at high enough Q 2 , can be viewed as the (elastic) scattering of the 
lepton from a quark or antiquark inside the proton mediated by the exchange of the 
corresponding virtual vector boson 7*, W* or Z*. Consequently, when the process 
is totally inclusive (one integrates over all final hadronic activity), can be fully 
described by two Lorentz invariants. With the momentum assignment of Figure 3 
one can construct the following invariants: 



s= {p + k ) 2 , t = (p-p') 2 , W 2 = (p + q) 2 , Q 2 = -{k' -k) 2 = -q 2 (10) 



defined as the centre of mass energy squared, the square of the four-momentum 
transfer between the proton and the final hadronic state, the invariant mass squared 
of the final hadronic state, and four momentum transfer squared between the lep- 
ton and the proton, respectively. It is also convenient to introduce dimensionless 
invariants (scaling variables) 



p_q 
p ■ k 



Q 2 



(ii) 



y = 



X = 



2p-q 



i.e. the inelasticity of the scattered lepton, and fraction of the proton momentum 
carried by the struck quark, or Bjorken variable, respectively. 
Notice finally that, ignoring masses, 



Q 2 = sxy W 2 = Q 



2 1-X 



(12) 



x 




r x(p') 



p(p) 



FIGURE 3. DIS kinematics. 



so that Q 2 max = s, the center of mass energy squared, and small values of x imply 
increasing W. As we have already said, in principle one is allowed to use any 
pair of these invariants to describe totally inclusive DIS. Usually x and Q 2 are the 
preferred ones. 

Another interesting point to remark concerns the resolving power of a DIS ex- 
periment. Clearly, the size d one can resolve inside the nucleon becomes smaller 
for large photon, or gauge boson in general, virtuality, namely 

he „ n GeV fm 

d ~ 0.2 — 

Q Q 

The 'magnifying power' is then d ~ 10~ 14 cm for Q 2 = 4GeV 2 , and d ~ 10~ 16 cm 
and for Q 2 = 40, 000 GeV 2 . 



Experimental Reconstruction 

Different DIS experiments cover different regions of the kinematical plane (x, Q 2 ) 
as shown in Figure 4 [9]. It is particularly interesting the very extended range 
covered by HERA: Q 2 ~ 0.2 - 10 4 GeV 2 and x ~ 10~ 5 - 10 _1 . 

The kinematics of the DIS events, namely the two invariants required to specify 
DIS processes, can be determined (particularly at HERA) from measurements on 
the electron (the lepton) alone, on the final hadronic state corresponding to the 
struck quark alone or on a mixture of both. In general, the preferred method 
depends on the particular kinematical region of interest, strongly correlated with 
the detector performance. 

The electron method: The input are the energy E' e of the final electron and the 
electron scattering angle 9 e measured with respect to the proton beam direction 
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FIGURE 4. Coverage of the (x, Q 2 ) plane by the various DIS experiments 



(9 e = 180° means zero electron scattering angle) 



y e = l-f sin 2 | Q 2 = 4£ e ^cos 2 | 

The hadron method: The input is the hadronic energy and the three-momentum 
components of the hadronic system. They are computed as the sum over all final 
state hadrons h 

_ Eh(E h -p z , h ) 2 _ V 2 x , h + pj h 
Vh ~ 2E e Qh ~ 1- Vh 

This method is also known with the name Jacquet-Blondel, because these authors 
were able to show that the contribution from hadrons lost in the beam pipe is 
negligible. 

The sigma method: This method is based on both electron and hadron measure- 
ments. The denominator of y^ is replaced by a sum over all final state particles 
including the scattered electron (certainly equal to 2 E e due to energy-momentum 
conservation). The name of the method comes from the introduction of the variable 
£ = T,h( E h - Pz,h) that appears in 

v = S Q 2 = E?sm 2 6 e ) 

^ E + ^(l-cos0 e ) ^ 1-yv 

Notice that the denominator in y^, is twice the true incident electron energy in the 
eventual case when the electron radiates photons that are not detected. 

There are other more sophisticated methods, well adapted to particular situa- 
tions, like the so-called double angle method or the PT method used by ZEUS 
[3]. 



Cross-sections 



The fundamental measurement in DIS experiments concerns the totally inclusive 
cross-section for 

£{k) + N(p) — > t'(k') + X(p') 

as a function of the kinematical variables defined above. For charged lepton-nucleon 
scattering mediated by the neutral current, the spin averaged cross-section is given 
in terms of the structure functions F 2 , F^ and F 3 



d 2 a± 4 nop' 



dxdQ 2 xQ 4 



(13) 



\- y +f) F 2 (x, Q 2 ) - £ F L (x, Q 2 ) T (y ~ tt! F 3 (x, Q 2 ) 



here a is the QED coupling constant. For Q 2 values below the Z° scale, the 
parity violating effects related to F 3 are negligible and all the process is due to 



7* exchange. Remember also that the longitudinal structure function F L , a QCD 
correction important for large y, is defined in terms of the standard F\ and F 2 , 
related to the transverse and longitudinal 7* N cross sections respectively, as 



F L (x,Q 2 ) = F 2 (x,Q 2 



1 + 



Q 2 , 



2xF 1 (x,Q 2 ) 



(14) 



and that in the naive quark-parton model, really valid at extremely high Q 2 , where 
quarks are considered free, massless, having spin 1/2 and without any p T developed, 
F L is zero because the Callan-Gross relation 

2xF 1 (x) = F 2 (x) 

is satisfied. Under this assumptions, the so called Bjorken scaling is fully satisfied, 
namely, structure functions are only function of the x variable. 
The ratio R is defined by 



R(x,Q 2 ) 



F L (x,Q 2 ) 



F 2 (x,Q 2 )-F L (x,Q 2 ) a T 



(15) 



which is obviously zero in the Callan-Gross limit and can be interpreted as the 
ratio of cross-sections for the absorption of transversely and longitudinally polarized 
virtual photons on nucleon. The differential cross-section (13) can be rewritten in 
terms of R(x, Q 2 ) as 



d 2 a 



Ana 2 



dxdQ 2 xQ 4 



Ml x 2 y 2 y 

1 -y „ n + ' 



1 + 



4M 2 N x 2 

Q 2 



Q 2 



2 l + R(x,Q 2 ) 



H^Q 2 ) 



(16) 



where the parity violating contribution was discarded. 

It is possible to specify the structure function F 2 in terms of the partons (or 
better quarks) within the nucleon as follows 



F 2 (x,Q 2 



x 



E 

/ 



q f (x,Q 2 )+q f (x,Q 2 ) 



(17) 



where the sum runs over all momentum distributions of quarks of the different 
flavours / excited (xqf(x, Q 2 )) and antiquarks (x cjf(x,Q 2 )) contained in the nu- 
cleon and e qf stands for the the different parton-photon couplings, namely their 
electric charges. Notice that the previous expression (17) is valid only in leading 
order of perturbation theory. In fact, there one has included the explicit Q 2 depen- 
dence which is the effect of implementing first order perturbative QCD. However, 
there is a particular renormalization and factorization scheme that can be used in 
higher order QCD, the so called DIS-scheme, where one can retain that expression 
at any order of perturbation theory. In other schemes, like the usual MS one, 



the expressions for the structure functions are more involved, including in general 
g(x,Q 2 ), the gluon density in the nucleon. 

Let us now refer for a moment to the inclusive reaction that goes via charged 
current 

Up + N^fT + X (18) 

where as usual we write N for an isoscalar nucleon. In the QCD improved parton 
model, the corresponding cross section reads 



d 2 a 2G 2 F M N E U ( M\ 



2 \ 2 
W \ 



dxdy 7i \Q 2 + M 2 V 



xq(x,Q 2 )+xq(x,Q 2 )(l-y 2 ) (19) 



M w stands for the charged intermediate boson mass. G F is the Fermi constant. 
The quark and antiquark distribution functions can be written in terms of the 
corresponding valence and sea flavour distributions in a proton as 

2 q(x, Q 2 ) = u v (x, Q 2 ) + d v (x, Q 2 ) + u s (x, Q 2 ) + d s (x, Q 2 ) (20) 

+2s s (x,Q 2 ) + 2b s (x,Q 2 ) 
2 q(x, Q 2 ) = u s (x, Q 2 ) + d s (x, Q 2 ) + 2 c s (x, Q 2 ) + 2 t s (x, Q 2 ) (21) 

One can introduce immediately the corresponding structure function F 2 for neu- 
trino scattering. It is of interest to compare this F% with the one for charged lepton 
DIS written above. Taking into account the sea corrections, the relation between 
both F 2 functions results 

F{ 18 V 5 q + q J ( ' 

This relation compares well with experimental data 



QCD evolution equations 

All the functions introduced above include an explicit Q 2 dependence. This fact 
reflects the presence of interactions among partons. Clearly, the parton model, 
that implies exact Bjorken scaling, has to be modified to include these interactions. 
Moreover, quarks and gluons are finally confined within hadrons by the non-Abelian 
gauge interaction in QCD. Even if QCD is with us for more than 25 years, it has 
not yet been possible to calculate in detail the structure of hadrons starting from 
quarks, gluons and their gauge interactions, making compatible the short distance 
almost free partons with the long distance confinement. This fact, related to the 
running character of the strong interaction strength that we have discussed above, 
is on the origin of our continuous use of structure (and fragmentation and frac- 
ture) functions in describing DIS. Fortunately enough, some factorization theorems, 



which are valid in QCD, allows one to split the problem into a perturbative part 
and a non-perturbative part. In other words, the cross-section can be expressed as 
the folding of the initial parton distribution function qf/N, a non-perturbative input 
given the density of partons / in the nucleon, with a lepton-parton cross-section 
that can be computed by means of perturbation theory inside QCD. This fact can 
be schematically written as 

% = ^2 Qf/N <S> o lq 
f 

and consequently, it can be sketched as in Figure 5 




N(p) 

FIGURE 5. Factorization in QCD 

It is clear that the non-perturbative part, namely, the part related to structure 
functions, has to be determined by fitting experimental data. Nevertheless, they 
are universal in the sense that once obtained from a given particular process, they 
can be used in connection with any other one. 

Regarding QCD as an improvement of the quark-parton model, we can say that 
the nucleon is not simply composed of three point-like quarks, the so called valence 
quarks but as soon as Q 2 increases, the vector boson (7*, Z* or W*) increasingly 
resolves the composition of the nucleon. One of the targets that the vector boson 
could find is one of the quarks, called sea quarks, which originate from a gluon, 
namely through 

9 -»• QQ 

being this gluon itself radiated from one of the valence quarks. In other words, 
increasing Q 2 means that the resolution of our "view" of the nucleon increases, so 
that the effective number of partons sharing the nucleon momentum also increases. 
Consequently, the probability of finding partons with small x has increased while 
the corresponding one to large x has decreased. This process clearly entails a 
violation of Bjorken scaling through a Q 2 dependence, that was born from QCD 
interactions. Consequently, to compare experimental data with QCD predictions, 
one has to compute perturbative QCD corrections to the fundamental process 



V*q -> q 



(23) 



Nowadays, these corrections, even if the calculation is very involved, are well known 
up to order a 2 [11]. We give here only an introductory discussion of the lowest order 
corrections. 

The relevant diagrams corresponding to the first terms of the perturbative ex- 
pansion in a s for the process (23) are presented in Figure 6 (We have omitted the 
vertex and self energy gluonic corrections, which are associated to the ultraviolet 
divergences. These are removed through coupling constant renormalization, and 
lead to the running coupling a s (Q 2 ) given in Equation (9)) 





FIGURE 6. Lowest QCD corrections to V* q -> q 

Because we are considering the inclusive case, the calculation implies phase space 
integrations which lead to divergent integrals when partons become collinear. They 
have to be regularized by means of dimensional regularization [12], a cut-off in 
the transverse momentum of the partons, or giving them masses. Any of these 
procedures introduce a regulator ji 2 which has the dimension of the energy scale 
Q 2 - 

For example, to order a s , quarks contribute to the structure function F 2 through 
the following convolution integral 

F« uarks {x,Q 2 ) = 




where qf(y) is the 'bare' quark density. The delta function term comes from the 
lowest order diagram in Figure 6, and represents the corresponding contribution 



to F 2 as in Equation (25). The first term between the curly brackets includes the 
factor \n(Q 2 /fi 2 ) which diverges as the regulator /x 2 goes to zero and is associated 
to the collinear gluon emission diagrams, while the second term collects the finite 
contributions from those diagrams. P qq (x/y) and R{x/y) are known calculable 
functions. 

The collinear singularities of course threaten the validity of the parton model, 
however they can be consistently removed from the partonic subprocess absorbing 
them into the 'bare' parton densities qf(y) defining renormalized parton densities 



(25) 



where M 2 is factorization scale chosen to separate the short distance ('partonic') 
effects from the long distance ('hadronic') ones. For the parton model to make 
sense, the renormalized parton densities must be process independent, i.e. must be 
the same for DIS, Drell-Yan, and any other process. Fortunately, this proves to be 
the case to all orders of perturbation theory, and one ends up with an expression 
for the (physical, and thus finite) structure function F2(x,Q 2 ) in terms of a finite 
renormalized parton density qf(x, M 2 ) and an also finite partonic cross section 



F« uarks (x,Q 2 ) 



x 



£4 r% f { y ,M*) 



(26) 



At this point, it is customary to chose the factorization scale M 2 equal to the 
energy scale Q 2 , factorizing the scale dependence of the cross sections into the 
parton densities. 

The crucial observation here is that although perturbation theory can not make 
an absolute prediction for qf(x,Q 2 ), from Equation (25) it follows 

<ki(x, M 2 ) 



dlogQ 2 



= as r 1 dy 
2-k Jx y 



x 



q l {y,Q 1 )P qq \- 



y , 



(27) 



which means that QCD actually gives the Q 2 dependence of the parton distribu- 
tions. 

By means of a similar procedure with the gluon densities, one can deal with the 
divergences in the diagrams initiated by gluons. However, due to the possibility 
of gluons emitting quark pairs which then interact with the photon probe, the 
evolution the evolution of the quark and gluon densities is no longer independent, 
but coupled. Taking into account all the 0(a s ) contributions, we end up with the 
so called DGLAP equations (or Altarelli-Parisi equations for short) 



dqj(x,Q 2 ) 
dlnQ 2 

dg(x,Q 2 ) 
dlnQ 2 



2^ 
2ir 



idy 

y 

*dy 

y 



qi(y,Q 2 )P qq 

T,<n(yiQ 2 ) p i 



X 



+ g(y,Q 2 )P t 



99 



X 



99 



-j+g(y,Q 2 )P gg 




(28) 



Collecting all the 0(a s ) contributions to F 2 (x, Q 2 ), we finally arrive to the most 
usual expression for the structure function: 



F 2 (x,Q 2 



)C f {x)® [q f (x,Q 2 ) + q f (x,Q 2 



(29) 



+C g {x) ®g{x,Q 2 ) 



where ® indicates a convolution integral and the index g refers to gluons. 

Deep Inelastic Scattering (DIS) experiments and hadron-hadron interactions pro- 
vide information on the parton structure of the nucleon and constraints the dynam- 
ics of the quark-gluon interaction given by QCD. In Figure 4 the kinematic regions 
in x and Q 2 for cross section measurements in DIS ep scattering, v scattering and 
jets in pp collisions were presented. Data from the last years [2] have a clear influ- 
ence on parton distributions and have motivated the update of parton distribution 
function analysis. Thes most recent sets are CTEQ5 [13], GRV98 [14] and MRST 
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FIGURE 7. F2 as a function of Q 2 and the corresponding QCD fit 
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[15]. The data playing a fundamental paper in these updates are the more precise 
ZEUS and HI determinations of F 2 P including F£ harm - the NMC and CCFR final 
muon-nucleon and neutrino data; the E866 pp and pd lepton pair production asym- 
metry; the W charge rapidity asymmetry; the DO and CDF analysis of inclusive 
single jet production and the E706 direct photon production. It should be remarked 
that both approaches differ in the selection of data for the global analysis that are 
sensitive to gluons at large values of x. This means that a better understanding of 
the gluon behaviour at large x is in order. 

The F 2 structure function measurements at HERA cover a wide range of five 
orders of magnitude in both x and Q 2 values as shown in Figure 4. One important 
point established by these data refers to the rise of the structure function with 
decreasing x. Moreover, QCD fits including NLO (next to leading order corrections) 
performed by the experimental collaborations are in very good agreement with the 
data even at low values of Q 2 ~ 1 GeV. In Figure 7, F 2 is presented as a function 
of Q 2 for fixed values of x from HI, ZEUS and the previous experiments, together 
with the QCD fit [2]. 

Figure 8 presents de neutrino and antineutrino cross section as a function of y 
and different values of x [2] . Data comes from the CCFR Fermilab experiment and 
is presented for (E v ) = 150 GeV. A good agreement with previous data is observed 
and the NLO QCD analysis is satisfactory. 

We finish this paragraph with a brief comment about some QCD technical details. 



The use of the renormalization group result for a s in the AP approach, specifically 
the replacement of a s by a s (Q 2 ) in Equation (28), has an important consequence. 
For each order of perturbation in a s where the factorization procedure is performed, 
a whole series of perturbative contributions is efectively re-summed, in addition to 
the contributions comming form the corresponding diagrams shown in Figure 6. 
This kind of contributions are depicted in Figure 9, the so called ladder diagrams. 



It can be shown that the improved AP approach takes into account contributions 
coming from the sum of all ladder diagrams in which the transverse momenta along 
the sides of the ladder are strongly ordered, namely: Ql -C k\ x <C • • • <C h\ { <C 
■ ■ • <C Q 2 . This condition, of course implies also a strong ordering for the momenta 
of the emitted partons. 

Notice that the calculation of such a n-th order ladder diagram implies integra- 
tions over the internal momenta. This nested integration can be carried out because 
of the kr ordering imposed, giving rise to results proportional to powers of a s (Q 2 ) 



Now, it is clear that large logarithms in Q 2 compensate the small values of 
a s (Q 2 ), that as we have seen, decreases logarithmically. Consequently, all graph 
with rungs up to n — > oo would have to be summed up. The so called leading 
log approximation (LO), which is equivalent to follow the renormalization group 
improved AP approach up to order a s , collects all the contributions proportional 
to 




FIGURE 9. A typical ladder diagram 



and In {Q 2 /Q 2 ). 




while the next to leading log approximation (NLO), equivalent to the second order 



extension of the AP approach, also includes contributions proportional to 



<{Q 2 ) 




n— 1 



The AP approximation is expected to be valid for Q 2 values sufficiently large but 
also for x not too small, just to ensure that small values of x do not give rise to 
other large logarithms. In other words, 



a s {Q 2 ) In (£)«a s (Q 2 ) In (|j 



< 1 



In is worth remarking that the AP equations can be solved analytically when a 
strong ordering in x is also requiered. In this way one ends with the double leading 
log approximation, (DLL), where the large logarithmic terms are of the form [16] 



Jt) h/ 1 



X 



This approximation is expected to be valid for large Q 2 and small x. 

At small values of x, the parton content of the proton is gluon dominated, 
this case the AP DLL can be obtained with the result 



In 



xg(x,Q 2 ) xg (x,Ql)exp, 
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In 



ln(Q 2 /A 2 ) 



ln(l/V) 



(30) 



that shows that the gluon density increases faster than a power of ln(l/a:). 

In some cases of interest, x is small enough but Q 2 is not sufficiently large to 
be inside the DLL regime. Under these circumstances, the AP approximation is 
not more valid. The BFKL (Balitsky, Fadin, Kuraev, Lipatov) equation has been 
proposed to tackle the limit behaviour of large 1/x and Q 2 finite and fixed. In this 
scheme, the Xi variables in the ladder are stongly ordered, namely: x <C x\ <C 
■ • ■ C X; < ■ • • C x; while there is no order on kx imposed. This approach ends 
with the so called leading log approximation in ln(l/x). The region of validity 
being 
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Due to technical reasons, the BFKL equation is written in terms of the function 
f(x,k^), related to the usual gluon density g(x,Q 2 ) that dominates at very small 
x by 
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The standard form of the equation is 
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FIGURE 10. A e// and dF 2 /d\nQ 2 . 
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This approximate BFKL equation can be solved analytically for fixed a s . The 
solution behaves like 

f(x, fc£) oc 
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with 
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(33) 



for iV c = 3 and a s = 0.19. 

Consequently, the gluon density rises like a power ofl/x for decreasing x, namely 

x g(x, Q 2 ) oc 

clearly faster than the AP DLL prediction (30). It should be noticed that a running 
a s and higher order corrections decrease the value of A. 

Assuming that gluons dominate at low x one can expect an approximate be- 
haviour for the structure function given by 

F 2 (x,Q 2 ) \ Q2 =cx- x *ff \ Q 2 

Figure 10 (a) shows the values obtained for A e // as a function of Q 2 from fits to 
ZEUS and E665 data. Figure 10 (b) shows dF 2 /dlnQ 2 as a function of x [17]. In 
both figures QCD and Regge inspired fits have been included. The schematic map 
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FIGURE 11. The schematic map of the plane ln(Q 2 ), ln(l/x) 



of the plane ln(Q 2 ), ln(l/x) shown in Figure 11 intents to divide the regions where 
each of the two approaches, DGLAP and BFKL, are in order [9]. 

In the figure, a saturation region is indicated, where one waits the gluon density 
to be so high as to prevent its continuous growing. 



In the analysis of hadron structure, more information can be obtained if one goes 
one step beyond totally inclusive DIS, namely to semi-inclusive DIS, where one of 
the final hadrons is also measured 



In describing semi-inclusive processes, in addition to quark distributions, the so 
called fragmentation functions are necesary. These functions [18] describe, or better 
parametrize, a given parton decay into a final hadron. 

It is clear that, in principle, the best process to analyze fragmentation functions 
is the semi-inclusive e + e~ annihilation, when a single hadron is fully detected in 
the final state 



Semi-inclusive 



e + h-y e' + h' + x 



e + e~ -> h(P) + X 




(34) 



It is also convenient to introduce ( that measures the fraction of the parton 
momentum carried by the final hadron. The corresponding cross section reads 



FIGURE 12. Schematic view of semi-inclusive e+ e_ 
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Here a l p is the parton % production cross-section, that in the parton model results 



ai(z,Q 2 ) = e 2 5(l-~z) 



(36) 



if i — q, q and vanishes if i — g. The function D^^Q 2 ) is the fragmentation 
function or probability for a parton i to decay into a hadron h carrying the fraction 
( of the parton momentum. 

When QCD corrections are in order, the fragmentation functions acquire a Q 2 
dependence. The evolution of these functions with the scale 



is given by the AP-like integrodifferential equations 

+ dy 
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D h q dt(y, t) P qq (z/y) + 2 / D h g dt(y, t) P gq (z/y) (38) 



(39) 



The functions Pij are exactly the same as those appearing in the AP equations for 
the structure functions (28). The corresponding graphical interpretation in terms 
of the basic QCD processes is included in Figure 13. 

It is worth noticing that the properties of the splitting functions are enough to 
guarantee the momentum sum rule 




FIGURE 13. QCD processes providing the evolution of fragmentation functions 

E CzD h q {z^) = l (40) 

h J0 

and the analogous one for the gluon fragmentation function. 

To solve the AP equations it is necessary to know the fragmentation functions 
at some value Q^. Exactly as in the case of structure functions, they have to be 
obtained from experimental data. In general one proposes a given form for the 
initial D-function, in agreement with data, and then AP equations provide the 
evolution with Q 2 . 



Fracture Functions 



For those who are not acquainted with fracture functions, let us briefly summarize 
that the main idea behind fracture functions is the realization that the most familiar 
perturbative description for semi-inclusive processes, based on parton distributions 
and fragmentation functions is, at least, incomplete [19]. 

In the usual approach to semi-inclusive DIS, the corresponding cross section is 
expressed by a convolution between parton distributions and fragmentation func- 
tions accounting the process in which the struck parton hadronizes into a detected 
final state particle. 

fe - x £:V + %> v)2 v eUtix ' Q2) x dHz - q2) (4i) 

Obviously, this approach only takes into account hadrons produced in the current 
fragmentation region, and what is more, within this approximation, and in leading 
order, hadrons can only be produced in the backward direction. Going to higher 
orders one finds a breakdown of hard factorization, as there are collinear singulari- 
ties that can not be substracted in the parton distributions or in the fragmentation 
functions. 

All this means that there are additional contributions missing, which are mainly 
target fragmentation processes and are included in the so called fracture functions 
Mf , defined by 



1 + (1 -y) 2 

2y* 



£e 2 (l-:r)Mf(*,z,Q 2 ) 



(42) 



dxdydz x (p + k) 2 



These functions, that can be thought as the probabilities to find a parton of a given 
flavour in an already fragmentated target, straightforwardly solve the factorization 
problem and also allow a LO description of hadrons produced in the forward direc- 
tion. The only subtlety regarding them is that they obey slightly different evolution 
equations; their scale dependence not only depends on their shape at a given scale 
but also on that of ordinary structure functions and fragmentation functions, re- 
flecting the fact that current and target fragmentation are not truly independent 
of each other. This is usually refered to as 'non homogeneus' evolution 



As for structure functions in totally inclusive deep inelastic scattering, QCD does 
not predict the shape of fracture functions unless it is known at a given initial scale. 
This nonperturbative information have to be obtained from the experiment, and, 
eventually, can be parametrized finding inspiration in nonperturbative models, as 
it is the case for ordinary structure and fragmentation functions. 

Fracture functions have been succesfully applied to the description of leading 
baryon production [20], have been extended to spin dependent processes [21], and 
as we shall see, are relevant for the description of diffractive DIS. 



It may seem strange that the next step in our study of parton structure relies in 
a particle, the photon, that it is not a hadron and of course has no structure by 
itself, as we very well know. However, a high energy photon can 'develop' structure 
when interacting with another object. In fact, as the electromagnetic field couples 
to all particles carrying the electromagnetic current, a photon can fluctuate into 
particle-antiparticle virtual states, particularly into quark- ant iquark pairs or even 
more complex hadronic objects with the same quantum numbers. As long a the 
fluctuation time is longer than the interaction time we can talk about the structure 
of the the photon, and deal with it as a 'very special' hadron. Notice that at high 
energies, the fluctuation of a photon into a state of invariant mass M can persist for 
a time of the order of r ~ 2 E 1 /M until the virtual state materializes by a collision 
or annihilation with another system [22]. 

There are many reasons for this peculiar hadronic structure to be considered 
'matter' of study. First of all, it gives the opportunity to study how QCD works in 
a completely different scenario. Theoretical studies of photonic parton distributions 
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FIGURE 14. Photon Structure. 



of real, i.e. on shell, photons have a long history initiated by Witten's work 
This description have been well established through measurements in 77 collisions 
at e + e~ colliders (PETRA, PEP, LEP) however, extending the idea of photon 
structure to high virtualities as in ep processes like dijet production, a new insight 
has been gained as we link DIS, photoproduction, and 77 interactions [24]. 

The concept of photon structure functions for real and virtual photons can be 
defined and understood in close analogy to deep inelastic lepton nucleon scattering, 
via the subproccess 7*(Q 2 )7(-P 2 ) — > X in e + e~ — > e ± X as in Figure 14, where we 
denote the probed target photon with virtuality P 2 = —p 2 by 7(-P 2 ) and reserve 
7*(Q 2 ) for the highly virtual one. The relevant differential cross section can be 
expressed, as in the hadronic case, in terms of the usual scaling variables x and y 
as 



dV(e 7 (P 2 ) — >eX) 2vra 2 s 



dxdy Q 4 

-y 2 F2 (p2) (x,g 2 ) 



(i + (i-y) 2 )Fl (p2 \x,Q 2 ) (44) 



with F^[ \x,Q 2 ) denoting the photonic structure functions. The measured e + e~ 
cross section is obtained by convoluting Equation(44) with the photon flux for the 
target photon 7(-P 2 ) [25]. The range of photon virtualities explored in the above 
mentioned experiment is given by 

m 2 e y 2 /(l -y)< P 2 mm < P 2 < P 2 max < ^(1 - y){\ - cos6 max ) (45) 

where m e is the mass of the electron, s e7 is the square of the c.m.s. energy, y is the 
energy fraction taken by the photon, and 9 max is the maximum scattering angle 
of the electron in the c.m.s. frame. Pmin,max are further determined by detector 
specifications and experimental settings. 

It is worthwhile noticing that if one could neglect photon fluctuations into com- 
plex hadronic states, then the dependence of \x, Q 2 ) in both x and Q 2 would 




be fully predictable in perturbative QCD. The pointlike process 7*(Q 2 )7(P 2 ) — ► 

qq, corrected by gluon radiation effects, yields a definite prediction for F^^ (x, Q 2 ). 
However, as it was anticipated, this is only the so called 'direct' contribution to the 
total process, named in opposition to the 'resolved' one, where the virtual photon 
7*(Q 2 ) strikes the nonperturbative 'dressing' of the photon fluctuation (Fig. 15). 

QCD corrections to the 'direct' component not only take into account the split- 
ting of partons into partons as in the ordinary AP evolution equations, but also the 
possibility of a photon to split into a quark-antiquark pair. This implies an inho- 
mogeneous term in the evolution equations and leads to a logarithmic enhancement 
in F^i \x, Q 2 ) which means positive scaling violations for all values of x (Fig. 16) 
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FIGURE 16. F 2 7 (x, Q 2 ) as a function of Q 2 . 



The 'resolved' component obeys the same evolution equations typical of the 
lepton-hadron interactions. In a Vector Meson Dominance approach for the photon 
fluctuations, the resolved component is expected to vanish like (1/P 2 ) 2 in connec- 
tion to the vector meson propagator, however it is not clear up to which values of 
P 2 the nonperturbative contributions are relevant. In recent years several sets of 
parton distributions for real and virtual photons have been proposed [27,28]. For 
virtual photons, different approaches have been followed in LO y NLO global fits 
to the available LEP data (Figure 17) [26]. 
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FIGURE 17. F 2 7 (x,Q 2 ) as a function of x. 



LEP data constraint reasonably well the quark distributions in the photon, how- 
ever, in order to constrain the gluon densitiy it is much more suitable to analyze 
dijet production processes in e + p collisions as measured at HERA. As we shall see, 
in this context the picture of 'direct' and 'resolved' is also more clearly illustrated. 
In this kind of process the parton content of the proton is used to probe the iqi 
partonic structure of the photon as in a hadron-hadron collision. The hard scale 
invoked to allow such inspection and guarantee the validity of a perturbative treat- 
ment, is the large transverse momentum of the final state jets. Figure 18 depicts 
LO 'direct' and 'resolved' contributions to dijet production e + p collisions. As it can 
be seen, even at the lowest order, photonic gluons contribute to the cross section, 
at variance with what happens in e + e~ collisions. 

In leading order the differential cross section for two jet production in ep collisions 





P (a) P (b) 

FIGURE 18. Direct (a) and Resolved (b) Contributions to dijet production. 



takes a very simple form when written in terms of the fraction of the photon energy 
intervening in the hard process, x 1 , the fraction of the proton energy carried by 



the participating parton, x p , and that of the electron carried by the photon, z [29] 



da 



dx 1 dx p dz dpr dP 2 



f, /e (z,P 2 )q^x„Q 2 ,P 2 )q^x p ,Q 2 ) 



da 



(46) 



Here, px is the transverse momentum of the jets, P 2 is the photon virtuality, and Q 2 
is the relevant energy scale of the proccess, taken in this case equal to p\. da/dpr 
represents the hard parton-parton and parton-photon cross sections [30]. 

The functions q 1 (x 1 , Q 2 ,P 2 ) and q p (x p , Q 2 ) denote the parton distribution func- 
tions for the photon and the proton, respectively. The first one reduces to 5(1 — rr 7 ) 
-the probability for finding a photon in a photon- for direct contributions, i.e. those 
in which the photon participates as such in the hard process. f y / e (z, P 2 ) is the un- 
integrated Weizsacker- Williams distribution [31] 
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(47) 



which has been shown to be a very good approximation for the distribution of 
photons in the electron, provided the photon virtuality is much smaller than the 
relevant energy scale [27]. 

The differential cross section in Equation(46) can also be written in terms of the 
parton pseudorapidities rji and r]2, which are constrained by the experimental set- 
tings, and the electron and proton energies E e and E p [32]. Both pairs of variables 
are related to the energy fractions by 
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(48) 



Kinematical restrictions constrain x 7 to lay in the interval \p\j [x v zE e E v \ x 
x p in \p 2 T /(zE e E p x™ ax ), 1] and z in \p 2 T /(E e E p ), 1]. 
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FIGURE 19. a) Event distribution measured in dijet photo-production, b) Ratio of 'direct' and 
'resolved' photon cross sections as a function of the photon virtuality. 



In the lowest order of QCD the 'direct' contribution is characterized by x 7 = 1, 
which means that all the energy of the photon participates in the hard interaction. 
Allowing gluon corrections, 'direct' contributions may come from x 7 7^ 1, however 
one should expect them concentrated arround a peak at high x 1 . Conversely, 
'resolved' contributions should be more copious al low x 1 as for ordinary hadrons. 
In fact, this trend is shown by the measurements, as can be seen in Figure 19a 
[33]. The variable x° bs used there to approximate x 7 is the jet level equivalent of 
Equation (48) 
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where pj) and r\i refer to the transverse energies and pseudorapidities of the mea- 
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FIGURE 20. Dijet cross sections for Q 2 < O.OlGeV 2 as a function of : 



sured jets instead of those of the partons. Of course, the relative abundance of 
'direct' and 'resolved' events depends strongly on the virtuality of the photon as 
can be seen in Figure 19b in accordance to the expectation about the suppression 
of the 'resolved' component as the virtuality increases [33]. 

Finally, it is quite instructive to anlyze the relative abundance of quark and gluon 
initiated contributions as a function of x 1 in a typical e + p dijet cross section. In 
Figure 20 data obtained by HI is compared to theoretical expectations based on a 
Montecarlo analysis using GRV parton distribution for the photon [27] and for the 
proton [34,35]. Notice that at large x 7 the 'resolved' component is clearly dominated 
by quarks, however as x 7 decreases, the photonic gluons become dominant. 



COLOUR SINGLET (POMERON?) 

In this section we draw our attention to a much more unfamiliar object, at least 
for the younger generation of physicists, although for some of us it may be an old 
acquaintance. We are talking about the pomeron, or in more modern language, 
about the colour singlet object which is exchanged when particles undergo strong 
interactions while preserving its nature. The possibility of studying the partonic 
structure of this old friend has been made real by a very recent generation of 
experiments and has added an extra quota of excitement in perturbative QCD. 

Before we plunge ourselves into the partonic structure of the pomeron, it would 
be very helpful to first fetch some of the tools used in past to sail in those waters. 



Peripheral model 

Before the advent of QCD, the presence of structures in the differential cross- 
sections was mainly explained by means of dynamical exchange mechanisms like 
in the so called peripheral model. In this approach it is found that there is a clear 
enhancement near the forward direction whenever the crossed channel of the reac- 
tion has the quantum numbers of a known particle. Then, peripheralism proposes 
to express the scattering amplitude as a sum of contributions corresponding to the 
exchange of a particle. Namely 

F — F n + F v + F B + ... 

with 



where it is clear that for small, physical (negative), values of t, there are peaks in the 
corresponding cross-section in s-channel, proportional to \F\ 2 . The indices 7r, V, ... 
represent pions, vector bosons, etc., and their masses hierarchize each contribution. 
The explanation of peripheral events cannot be found in QCD, since the exchanged 



particle is not a gauge boson that can be treated perturbatively. In fact the particles 
involved are composite objects such mesons or baryons. 

When the exchanged particle has a spin J, the previous expression for the am- 
plitude has to be replaced by 

j-i Pj(cos9 t ) 
Fi oc — — 

mj — t 

and, if one remembers that the scattering angle in the t-channel reads 

2s 



cos 6+ — 1 + 



t — Am 2 



it is found that the numerator of Fj behaves asymptotically as s J for fixed t. This 
is not tenable because the Froissart bound is not satisfied for J > 1. Moreover, 
that expression is purely real and does not satisfy the principle of analyticity. 

There is a way out of the above mentioned difficulties, while maintaining the 
spirit of an exchange mechanism, now of an entire family of related particles. This 
proposal is based on the Regge poles ideas. These ideas emerge naturally in po- 
tential scattering theory through the analyticity properties of the amplitude. The 
Regge pole model [36] is based then, upon the assumption of the following behaviour 
for the scattering amplitude 

lim F(s, t)~Y PiCt) , , X1 I - ) (50) 
♦oo, fixedt K,J Y sin[7rai(t)] Vs / 
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here & = ±1 is called the signature of the corresponding Regge pole. The functions 
oti{t) are called Regge trajectories and /3i(t) residues of the poles that in general 
factorize between the initial and the final channels. Both functions ««(£) and /3j(i) 
are analytic functions of t. 

The high energy model so introduced is understood as an exchange model where 
the Regge trajectories a^t) are such that they pass through the spin values of the 
particles, or better the resonances, each time positive t takes the values of their 
mass squared. In this way, the collective effect of the exchange of all members of 
the family is taken into account. Clearly, the Froissart bound (total cross-sections 
cannot increase faster than In 2 s) is satisfied as soon as 

ai(t < 0) < 1 

in the space-like region of negative t, relevant for the sacttering process. 
Let us take as an example the case of pion-nucleon charge-exchange 

7T~ p — > 7T° Tl 

that has always been the paradigmatic case of reggeology. The t-channel of this 
reaction is 

7T~ 7T° — > p U 



having the quantum numbers; Q = 1, I = 1, Y = 0. G — +1 that correspond 
to the well known meson resonances p(770MeV) of spin 1 and #(1680 Me^/) of 
spin 3. This sequence of resonances defines precisely the meson Regge trajectory of 
isospin 1, G-parity +1 and positive signature that is noted a p . The continuation of 
this trajectory to negative t values, takes us into the physical region corresponding 
to the s-channel charge-exchange reaction. The parameters of a standard linear 
trajectory are the intercept at t — and the slope. In the present case, as it is 
shown in Figure 21, they are 




FIGURE 21. p-Regge trajectory and it p — > 7r° n data 



a p (0)~0.55 a' ~ 0.9 GeV' 2 
respectively. Consequently, the differential cross section has the behaviour 

~ s 2 (Mo)-i) „ s -i 

t=0 

But what happens if the crossed channel of the reaction, and correspondingly the 
eventual exchanged object, has the quantum numbers of the vacuum? The answer 
arrives in the following paragraphs. 

Diffraction 

Landau and his collaborators, in the fifties, introduced the term diffraction in 
high energy physics used in complete analogy with the well known phenomenon in 
Optics that occurs when light interacts with obstacles or holes whose dimensions 
are of the order of the electromagnetic radiation wavelength. The interaction of a 
hadron could be think as the absorption of its wave function caused by the different 
channels open at high energy and consequently the name diffraction seems to be 
intuitive. For a recent review of this topic see [37]. 

In the Fraunhofer limit, namely when the product of the wave number times 
the area £ of the obstacle is of the order of the observation distance, the energy 




distribution at the observation point, given by the Kirchoff formula, can be written 

as 

T(x,y,z)^— — / d 2 bS(b)e^ b 

where b stands for the impact parameter, r is the position of E, \q\ = k sin6 is the 
2-D momentum transfer and the scattering matrix is expressed as 

s(b) = i - r(b) 

in terms of the so called profile function T(b). From the expression for T, one can 
immediately obtain the scattering amplitude 

M = ¥- [d 2 br(b)e^ b 

i.e.: given as the 2-D Fourier transform of the profile function. 

When the function T is spherically symmetric, the integral becomes a Bessel 
transform, namely 

f{q) = ik bdbT(b)J (qb) 
Jo 

meaning that if the profile function is merely a disk of radius R, the amplitude 
reduces to the black-disk form 

qn 

The prediction coming from this very simple model is a series of diffractive maxima 
and minima, entirely similar to the case in Optics, and have been clearly observed 
in several experiments 

In the specific field of particle physics, diffraction is said to be the dominant 
process of scattering at high energy if no quantum numbers are exchanged between 
the colliding particles. In other words, diffraction dominates asymptotically as 
soon as the particles in the final state have the same quantum numbers of the 
incident ones. This sort of definition of diffraction clearly includes, for the two body 
scattering, three cases: elastic scattering, single diffraction and double diffraction. 
In the first process the outgoing particles are exactly the same as the incident ones. 
In the second case, one incident particles goes out unmodified while the second 
one gives rise to a resonance, or to a bunch of final particles, with total quantum 
numbers coincident with its own ones. Finally, when double diffraction occurs, each 
incident particle gives rise to a resonance, or to a bunch of final particles, with the 
same quantum numbers of the initial ones. 

Pomeron 

It is possible to discuss diffraction from the viewpoint of an exchange model, in 
particular within the framework of reggeology. Remember that we have found 



(51) 



precisely that an enhancement in the differential cross section can be expected 
whenever the quantum numbers of the crossed channel correspond to an existing 
particle. However, in the elastic case at high energies, all the hadronic systems show 
the diffraction peak while the exchange of a given particle gives rise to different 
contributions in different systems. Moreover, via the optical theorem, the forward 
elastic amplitude is connected to the total cross section and this fact cannot be 
understood in terms of only one exchange process. Nevertheless, the diffraction 
peak can be interpreted in terms of a special trajectory that summarizes all the 
diffractive contributions and has the vacuum quantum numbers: the pomeron. Its 
contribution to the amplitude is 

1 _L e -ina F (t) / „ \ a (t) 

F dif (s,t) = -p„(t) (-) (52) 

sm[7rci; 2P (t)J \s J 

Consequently, the asymptotic behaviour of the total cross section (via optical the- 
orem) is given by 

so 

if the intercept of the pomeron trajectory is the maximum value allowed by the 
Froissart bound, namely ^(0) = 1. This fact implies that the total cross section 
in this approximation behaves asymptotically as a constant. It is clear that log- 
arithmic corrections are always allowed. Moreover, they are necessary in order to 
cope with experimental data. 

The pomeron concept answer the question at the end of the previous section 
using the reggeon jargon: a diffraction process is dominated by the exchange of a 
pomeron, that in fact means exchange of no quantum numbers. 

It is interesting to mention that Donnachie and Landshoff [38] were able to de- 
scribe all available total cross section data on pp, pp, up and Kp by a very simple 
Regge-Pomeron inspired parametization of the form 

atot =Xs 0.0808 +rs -0.4525 

where X and Y are reaction dependent parameters. Clearly the first exponent of 
the center of mass energy squared s correspond to a pomeron intercept of a P (0) = 
1.0808, while the second exponent comes from a typical Regge intercept 0:^(0) = 
0.5475. The successful fit is shown in Figure 22. 

Trying to connect the pomeron with the language and understanding of QCD, 
one can imagine it as a colour singlet combination of partons such as the simplest 
picture of a pair of gluons proposed by Low and Nussinov [39] . Two-gluon exchange 
is compatible with all soft phenomenology except that this simple model gives a 
constant, not a rising, total cross-section. Clearly in order to analyze the parton 
content of pomeron, a hard or high-momentum transfer interaction analogous to 
the one used to find the quark-gluon content of the proton is necessary. This 
kind of processes is called hard diffraction making reference to the hard scale that 
ultimately allows the pertubative description. 




Nowadays, people is trying to go beyond in the understanding of pomeron prop- 
erties following different theoretical approaches [40], although there remain several 
open questions. The main point to be solved is related to the precise relation be- 
tween hard and soft diffraction. As we said, to study the pomeron parton content, 
small distances and high-momentum transfer are necessary, while the natural envi- 
ronment of the soft pomeron is at large distances and low-momentum transfer. In 
other words, the clear notion of a pomeron is still missing. 

Most of the recent renewal of interest in diffraction was triggered mainly by a 
theoretical proposal of Bjorken [41], who pointed out a new signature of diffraction 
related to the presence of large rapidity gaps. These rapidity gaps were found both 
at HERA and at the Tevatron. 

The presence of diffractive processes at HERA should be expected since the 
photon behaves like a hadron in several circumstances. In fact, about 40 % of the 
photoproduction events are of diffractive character. In this kind of processes, the 
pomeron is exchanged and the incident proton remains as a proton or is diffractively 
dissociated in a state with the same quantum numbers. Consequently, as it will 
be discussed below, there appears a large rapidity gap between the proton, or the 
diffractive system, and the hadrons coming from the system into which the photon 
was diffracted. The quite unexpected fact observed at HERA was the presence of 
large rapidity gap events also in the DIS domain [42]. The observation that also a 
highly virtual photon can participate in a diffractive process, opened the road to 
the analysis of the pomeron structure. In fact, if Q 2 , the virtuality of the photon, 
is larger than say 4 GeV, the scattering 7*-pomeron can be treated perturbatively. 
The impact of this DIS diffractive data is also evident because one is facing in a 
DIS process characterized by a large Q 2 scale, diffractive properties which were 
been expected at a soft scale. 



Large rapidity gaps 

As it was predicted by Bjorken [41], the most evident signal for diffractive physics 
at high energy, is the presence of large rapidity gaps. Just to get some insight on 
this jargon, remember that a single inclusive diffractive reaction, noted as 

A(k A ) + B(k B ) ^ A'(k' A ) + X(k x ) 

implies no exchange of quantum numbers different from those of vacuum. In this 
inclusive case, at variance with the elastic reaction, besides the two variables (k 
and 9 or s and t) a third variable is needed to describe the process. Generally it is 
used 
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or, alternatively, the Feynman xp defined as 
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Another useful variable is rapidity, defined by 

E'a + ^A long 



y = In 



E'a ~ k Along 



(55) 



or equivalently, the pseudorapidity which is its limit for large energy 

T) = — In (tan — J 



(56) 



In a typical DIS event, the struck parton of the proton, emerges forming an angle 
9 a' with the proton remnant direction, as it is shown in Figure 23. 




FIGURE 23. DIS event diagramm 



This angle can be expressed in terms of the difference in total pseudorapidity r\ 
between these directions, namely 




on 



If the center of mass energy squared of the system is s, the pseudarapidity interval 
covered is 



where the first term correspond to the pseudorapidity covered by the 7* — p system 
and the second to that covered by the 7* — quark system respectively Clearly x is 
the Bjorken variable that measures the amount of momentum of the proton carried 
by the quark. 

The confinement property of QCD can be rephrased by saying that the struck 
quark and the proton remnant are connected via a colour string in order to end 
with colourless hadrons. For this reason, the mentioned pseudorapidity gap A77 is 
filled with particles during the hadronization process. Moreover, when the value of x 
decreases, the average hadron multiplicity increases, making less likely the visibility 
of any rapidity gap in the DIS event. In principle, the rapidity gap between the 
remnant proton, and the struck quark jet is then exponentially suppressed. 

After this kinematical preface, we briefly discuss the so called leading particle 
effect, characteristic of high energy hadronic interactions. Around 10% of the 
inclusive hadron scattering events present a Lab system configuration where the 
incident particle flies apart essentially unscattered in almost the forward direction 
leaving behind a stream of slow moving produced particles. The experimentally 
determined cross section looks entirely similar to the elastic case. Namely, between 
the initial particle and the final fast one, there is no change of quantum numbers 
and the reaction is clearly of diffractive type. There is only a minor loss of mo- 
mentum to produce the slow particles. This effect requires the fast particle to be 
exactly the same as the incident one. 

When the leading particle, the one we called above A', is produced diffract ively, 
and because the process has this characteristic, one has to expect the remaining 
cluster X to be produced at the opposite end of the rapidity spectrum. That is 
usually expressed by saying that the reaction presents a large rapidity gap from 
A'. Diffraction definition makes evident this concept. In fact, as A' has exactly 
the same quantum numbers as A, no quantum numbers at all can be exchanged 
until X is produced. If a particle were produced in between, it would mean that 
some quantum number has been exchanged and the process is no longer diffractive. 
Consequently, in this large rapidity gaps events, as no particle is produced between 
A' and X, they are clearly separated in rapidity. 




Consequently, 




Large rapidity gap events can be generated in soft processes as elastic and single 
diffractive hadron-hadron scattering, where the momentum transfer t between the 
scattered particles is small, and the center of mass energy s large, i.e., peripheral 
scattering. 

It should be stressed that nowadays, diffraction is considered synonym of large 
rapidity gap processes. In other words, the evidence of diffractive hadronic events 
comes mainly from large rapidity gaps as have been observed in DO and CDF 
Collaborations at the Tevatron [43]. 

On the other hand, at HERA, the electron-proton collider, diffraction has also 
been extensively observed using both the rapidity gap criteria for diffraction and 
also the detection of final state protons. In semi-inclusive DIS, whenever the de- 
tected final hadron coincides with the incident one, one is dealing with diffraction. 
In order to follow the outgoing proton a Leading Proton Spectrometer (LPS) was 
designed. 

Hard Diffraction at HERA 

Before discussing the interesting HERA DIS diffractive events, let us consider the 
particular semi-inclusive processes where the detected final hadron exactly coincides 
with the initial one, namely 

£(k e ) + N(k N ) -> £(k' e ) + N(k' N ) + X(k x ) 

In this case the scattering diagram looks like the one in Figure 24 
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FIGURE 24. Diffractive DIS event diagramm 

The process is clearly diffractive since no quantum numbers are exchanged be- 
tween the 7* and the nucleon N. As a consequence, one has for the remnant X 
the quantum numbers Jj^ c = 1 exactly equal to those of the incoming photon. 
Notice by the way that the production of any vector boson, the explicit replacement 
of X by any V in the reaction above, is a diffractive process because the V meson 
quantum numbers are precisely those of 7*. 



The kinematics related to Figure 24 needs further variables to be introduced. 
They are usually defined as 

t=(k N - k' N ) 2 , M 2 X = (k N — k' N + qf , My = k'x (57) 
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Notice that frequently £ is used instead of x p and that xf = 1 — x F . Clearly, the 
range of values of X j X jp . X J? and ft is (0, 1). Moreover, ft could be understood as the 
momentum fraction carried by the parton directly coupled to 7*. 

In the context of diffractive DIS, it is customary to define diffractive structure 
functions, in analogy with the ordinary ones, through the expression of the corre- 
sponding cross section 
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where 



if (4) = F^ (4) + F L D(4) 



It is also of interest the t-integrated diffractive structure function defined en each 
case by 



F^\x F ,(3,Q 2 )= f Wm "d\t\F^\x p ,p,^,t) 

** 5 m.in. 



(60) 



with % — 2, T, L. The superindices 3 and 4 refer to the number of variables in the 
structure functions. In the integral \t\ min is the lower kinematic limit of \t\ and 
\t\max has to be specified in each case. 

In the most naive Regge inspired approach, the diffractive structure function 
is assumed to be given by the product of the probability f p / p (x p ,t) to find a 
pomeron in the incoming proton, which only depends on the variables x p and t, 
and a pomeron structure function '((3,Q 2 ), given by parton densities which are 
assumed to behave according to Altarelli-Parisi evolution equations and factorize 
as ordinary parton distributions [44]. 

In recent years different theoretical approaches have been proposed for the de- 
scription of hard diffraction and specifically for the diffractive structure functions. 
See for example [40] and references therein. 

Some of these approaches modelize the diffractive interaction in the proton rest 
frame as an exchange of two gluons between the proton and the hadronic system 
into which the photon has evolved as depicted by the upper line diagramas in 
Figure 25. Others, also in the proton rest frame, evaluate within a semiclassical 
analysis the effect of the colour field of the proton into the photon hadronic system 
as in the lower line diagrams of Figure 25. The main difference between the above 
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FIGURE 25. Different approaches for diffractivc DIS 



mentioned models is that in the former only two gluons are exchanged while in the 
latter a multiple exchange of soft gluons is assumed [40] . 

A third strategy consists in dealing with the diffractive interaction as a spe- 
cial kinematical limit of a more general semi-inclusive process, i.e. regarding the 
diffractive structure function just as the low x p limit of the fracture function of 
protons into protons [19]. Doing this, the whole perturbative techniques can be 
rigurously applied without need to make additional assumptions in a program sim- 
ilar to what has been done for structure and fragmentation functions [45]. Another 
advantage of this last strategy is that in this framework the large x p behaviour of 
the diffractive structure function can be explored with leading proton production 
experiments. The fracture function approach can be in some sense related to the 
two gluon exchange model by a frame transformation as it is sketched in Figure 26. 

A remarkable thing to notice regarding these kinds of models for diffraction and 
also the factorization approach, is that although they seem to differ in rather strong 
assumptions, they all give a reasonably good account of the data suggesting a large 
gluon content in the diffractive structure function with gluons concentrated at high 
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FIGURE 26. Correspondence between the two gluon exchange models and the factorization 
approach. 
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FIGURE 27. Diffractive DIS data. 



In Figure 27 measurements of the diffractive structure function F 2 (x p , (3 , Q 2 ) 



a function of x T 



obtained by HI at HERA are compared with the fit coming 
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FIGURE 28. Scale dependence of F 2 {x IP ,f3, Q 2 ). 



from a factorization approach [45]. As for ordinary structure functions, QCD pre- 
dicts the scale dependence of the diffractive structure function. Figure 28 shows 
the agreement between the data and the behaviour expected in the factorization 
approach. 

Having established a rigorous and precise description of diffractive DIS we can go 
back to the most naive Regge inspired approach and see which of their hypothesis 
or assumptions may survive. The first thing to notice is that although in prin- 
ciple diffractive parton distributions or fracture functions obey non homogeneus 
AP evolution equations, in the low x p limit the non homogeneus contributions 
are numerically negligible so the standard assumption is a very good approxima- 
tion. Instead, what fails, at least in the unrestricted kinematical range accesed by 
HERA, is the flux factorization hypothesis, i.e. the posibility to factorize the x p - 
dependence of the structure function as a simple power of x F . Global analysis of 
the data shows not only deviations form this behaviour, suggesting the admixture 
of Regge exchanges, but also a /9-dependent level of admixture [45]. 

HERA also measures other diffractive procceses like diffractive production of vec- 
tor mesons and diffractive photoproduction of jets [46]. This last kind of processes, 
depicted in Figure 29, is particularly exciting as it represents a combined test of 
our knowledge of the parton structure of the photon and the colour singlet. ^From 
the theoretical point of view it is also very important because the resolved photon 
contribution to the process is in fact a hadron-hadron diffractive process for which 
hard QCD factorization may be broken. 

Indeed, although hard factorization can be proven to be valid for diffractive DIS, 
the proof fails for hadron-hadron interactions. The posibility of soft interactions 
taking place before the hard scattering occurs may spoil factorization in hadron- 
hadron collisions [48] , although there are no model-independent estimates of how 
large the factorization breaking effects may be or under which circumstances can 
be found. 

Hard factorization has been a very active theoretical topic in the three preceding 
years, and it is particularly relevant for the discussion of Tevatron diffractive data. 
Figure 30 shows a comparison between HI dijet photoproduction data and the 
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FIGURE 29. Direct and Resolved Diffractive Dijet Photoproduction 



corresponding prediction computed with photonic and diffractive parton distribu- 
tions [47]. The comparison suggests that for diffractive photoproduction the fac- 
torization breaking mechanisms are either neglible or beyond the accuracy of the 
present data. 
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FIGURE 30. HI dijet photoproduction data and the prediction computed with photonic and 
diffractive parton distributions 




Hard Diffraction at Tevatron 

As we mentioned earlier, large rapidity gaps in proton-antiproton collisions have 
been observed by both DO and CDF collaborations at Tevatron. The typical 
rapidity gap events observed at Tevatron are classified acording to three categories: 
hard single diffraction, hard double pomeron exchange, and hard colour singlet. 

In the first case, Figure 31a, a large rapidity gap in the forward direction is found 
between the outgoing antiproton and the debris produced by the proton-pomeron 
interaction. In the double pomeron exchange, Figure 31b, both the proton and the 
antiproton survive the interaction and emerge leaving rapidity gaps between each 
of them and the debris of the pomeron-pomeron interaction. In the third case, both 
the proton and the antiproton disociate but their corresponding remnants leave a 
large rapidity gap between them, Figure 31c. 

For a straightforward comparison with HERA data the most simple topology to 
analyze is hard single diffraction, as the only change to be made in the conceptual 
framework developed in the previous section is the replacement of the lepton probe 
(the positron) used at HERA by an hadronic one (the proton) that takes its place 
at the Tevatron. The diffractive structure function for the antiproton is simply 
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FIGURE 31. Rapidity gap topologies at Tevatron. 



related by charge conjugation to that for the proton, and which is measured at 
HERA. Whithin this topology, several observables can be measured such as dijet, 
W ± , and Z° production. 

Performing this kind of analysis what has been found is that even though the 
Tevatron data allow a description in terms of diffractive parton densities, the pre- 
dictions computed with those coming from HERA largely overshot Tevatron data 
signaling a breakdown of factorization [49], which remains to be fully understood. 

As we have already pointed out in the previous section, the factorization breaking 
seen Tevatron is not at all unexpected, since QCD factorization may be broken in 
hadron-hadron diffraction due to soft exchanges. In addition to these, there are also 
some other effects which conspire against factorization and need some consideration. 

For the analysis of these effects it is crucial to take into account that the standard 
criterium used to identify diffractive events both at HERA and at Tevatron, and 
thus the correspoinding contributions to the diffractive parton distributions and 
the Tevatron observables, respectively, is the presence of rapidity gaps of a given 
size. Obviously, changing the size of these rapidity gaps, the number of events 
selected also changes, effect that in fact has been observed and studied by the 
different experiments. In Figure 32 we show as an example dijet photoproduction 
events generated with a Monte Carlo in the kinematical regimes of HI and ZEUS 
experiments as a function of the pseudo-rapidity rfj^ of the most forward particle 
belonging to the system X [47] . Events to the left of the thick vertical line (rapidity 
gap limit) are those taken into account while those to the right are discarted. 

It is worthwhile noticing that even for the same kind of process, the different 
kinematical regimes covered by both experiments given for example by cuts in the 
transverse energy E T and pseudo-rapidity r\ of the two most energetic jets, etc., 
yield rather different distributions in f]^ x and thus, are affected in a different way 
even for the same pseudo-rapidity gap definition. For ZEUS, approximately 25 % of 
the total number of events are concentrated at rj^ x < 1.8, while for HI kinematics 
approximately 10 % of them survives the same constraint. 

Clearly, the situation is much more involved if we try to relate pseudo-rapidity 
gap data from different processes as we are doing in the HERA- Tevatron compar- 
ison., In other words, is not evident what is the relation between DIS events with 
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FIGURE 32. Dijct photo-production events generated in the kinematical regimes of the a) HI 
and b) ZEUS experiments 

a pseudo-rapidity gap of given size, and proton-antiproton events with pseudo- 
rapidity gaps of a different size, and obtained in a completly different kinematical 
regime. 

It is somewhat paradoxical that the same rapidity gaps that helped so much the 
development of hard diffraction, seem at this point to be a burden. In any case, 
ongoing experimental programs based on the Roman Pot technique, which allows 
a positive identification of the final state hadron as signature of the colour singlet 
exchange, will surely allow a much more deep and complete picture of these exciting 
phenomena. 

CONNECTIONS 

Throughout the present lectures we have review the most significant features of 
partons, the fundamental protagonists of QCD, as seen in the three rather different 
enviroments, the proton, the photon, and the pomeron (colour singlet exchange). 
As we have seen, and in fact was anticipated in the introduction, these objects, 
or more precisely the high energy processes involving them, are three of the main 
benchmarks of perturbative QCD. 

In the process of developing the predictive power of QCD we have introduced 
structure, fragmentation, and fracture functions, objects that link intimately both 
our knowledge and our ignorance about the partonic structure. 

The fundamental prediction of QCD in these enviroments, the energy scale de- 
pendence of the corresponding cross sections or structure functions, is beautifully 
illustrated in Figure 33 [50]. The observed behaviour in each case can be traced 
back to a basic parton interaction mechanism which alternatively dominates over 
the others: gluon radiation, quark pair creation from gluons, and from photons, 
respectively, and to the respective parton compositions. 

At intermediate x, the case in the figure, the proton structure function is domi- 
nated by valence quarks, however increasing the resolving power of our probe, gluon 




FIGURE 33. Q 2 Dependence of the proton, the color singlet, and the photon. 

radiation depletes dramatically these densities through the corresponding term in 
the AP evolution equation. At variance with this scenario, the diffractive structure 
function is dominated by gluons at intermediate /3, and consequently our lepton 
probes mainly quarks produced by pair creation from gluons, and the corresponding 
evolution. The logarithmic increase in the photon structure function is caused by 
quark pair creation from photons, which in the lowest order approximation appears 
in the AP equations as Q 2 independent term. 

Finally, the partonic structure as seen in our benchmarks in some way relate 
the physics made in the three main high energy physics laboratories. LEP pro- 
vides us with the most precise determinations of photon structure functions and 
also fragmentation functions, using electron-positron collisions. Tevatron in turn, 
contribute to our understanding of the proton structure and also to diffraction 
by means proton-antiproton interactions. Positron-proton collisions at Hera, not 
only test the structure of the proton and diffraction by themselves, but also brings 
toghether the proton and the photon structure functions in inclusive photoproduc- 
tion, relates the latter with fracture functions in the diffractive photoproduction, 
and explores the relation between fragmentation and fracture functions in leading 
baryon production. 
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